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Spin— orbit resonance, transit duration variation and 
possible secular perturbations in KOI-13 



o 

(N 
C 

0\ 



Gy. 



M. Szabo 1 ' 2 *, A. Pal : ' 3 , A. Derekas 1 , A. E. Simon 1 ' 2 , T. Szalai 4 , L. L. Kiss 1 ' 5 

1 Konkoly Observatory of the Hungarian Academy of Sciences, PO. Box 67, H-1525 Budapest, Hungary 
2 Department of Experimental Physics and Astronomical Observatory, University of Szeged, H-6720 Szeged, Hungary 
3 Department of Astronomy, Eotvos University, Pdzmdny Peter setdny 1/A, 1117 Budapest, Hungary 
^Department of Optics and Quantum Electronics, University of Szeged, Dom ter 9., H-6720 Szeged, Hungary 
5 Sydney Institute for Astronomy, School of Physics A28, University of Sydney, NSW 2006, Australia 



Accepted Received; in original form 



in 
6 



> 

m 

^1" 



ABSTRACT 

KOI-13 is the first known transiting system exhibiting light curve distortions due 
to gravity darkening of the rapidly rotating host star. In this paper we analyse publicly 
available Kepler Q2-Q3 short-cadence observations, revealing a continuous light vari- 
ation with a period of P ro t = 25.43 ± 0.05 hour and a half-amplitude of 21 ppm, which 
is linked to stellar rotation. This period is in exact 5:3 resonance with the orbit of KOI- 
13.01, which is the first detection of a spin-orbit resonance in a host of a substcllar 
companion. The stellar rotation leads to stellar oblateness, which is expected to cause 
secular variations in the orbital elements. We indeed detect the gradual increment of 
the transit duration with a rate of (1.14 ± 0.30) x 10 -6 day/cycle. The confidence of 
this trend is 3.85-cr, the two-sided false alarm probability is 0.012%. We suggest that 
the reason for this variation is the expected change of the impact parameter, with a 
rate of db/dt = —0.016 ± 0.004/yr. Assuming b « 0.25, KOI-13. 01 may become a non- 
transiting object in 75 — 100 years. The observed rate is compatible with the expected 
secular perturbations due to the stellar oblateness yielded by the fast rotation. 
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1 INTRODUCTION 

KOI-13 (KIC 009941662) is a unique astro-physical labora- 
tory of close-in companions in an oblique orbital geometry. 
The system consists of a widely separated common proper 
motion binary of A-type stars, one hosting a highly irra- 
diated planet candidate with P or b ~ 1.7626 day (Borucki 
2011). The host star of KOI-13.01 is the brighter compo- 
nent, KOI-13 A, and rotates rapidly (t>sini as 65 — 70km/s, 
Szabo et al. 2011). The transit curves show significant dis- 
tortion that is stable in shape, and it is consistent with a 
companion orbiting a rapidly rotating star exhibiting grav- 
ity darkening by rotation (Barnes 2009). Barnes et al. 2011 
derived a projected alignment of A = 23° ± 4° and the star's 
north pole is tilted away from the observer by ip = 48° ± 4° , 
therefore the stellar inclination, i* = 43° ± 4° (assuming 
M* = 2.05 Mq). The mutual inclination is <p = 54-56°. 
Companion is determined to have a mass of 9.2 ± 1.1 Mj 
(Shporer et al. 2011) and between 4 ± 2-6 ± 3 Mj (Mazeh 
et al. 2011), and a radius of 1.44 Rj (Barnes et al. 2011). 
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On 23 September, 2011, new Kepler photometry (Short Ca- 
dence Q3 data) became available. Here we suggest that the 
photometric data reveal the stellar rotation, and give ob- 
servational evidence for transit duration variations (TDV) 
which are a sign of secular perturbations. 



2 ROTATION OF THE HOST STAR 

After removing the orbital phase-dependent light variations 
from the Kepler photometry of KOI-13, we detected an addi- 
tional periodic signal that we argue below is due to rotation 
of the host star. The light curve was corrected with a di- 
lution factor of 1.818 (Szabo et al. 2011). The folded light 
curve with P or b is plotted in the upper panel of Fig. [JJ and 
we see the ellipsoidal and beaming variations (Mazeh et al. 
2011, Shporer et al. 2011). The reference flux value is taken 
to be the centre of the secondary eclipse, when only KOI-13 
A is visible. We smoothed and subtracted this variation from 
the out-of-transit light curve, and de-projected the residuals 
into the time domain (upper middle panel). A prominent 
frequency was detected in these residuals at 0.9437 1/day, 
i.e. at 25.43 ±0.05 hour period (Fig. [TJ lower middle panel). 
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Figure 1. Top: out-of-transit ligh tcurve of KOI-13, corrected 
for the third light. The reference flux level is KOI-13 A, i.e. the 
flux measured at eclipse. After subtracting this signal from the 
time series (middle upper panel), a signal with 25.432 hour period 
is detected with high significance (lower middle panel). Bottom: 
The phase diagram with 25.432 h period, after removing the orbit- 
dependent variations. 
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Figure 2. Fourier-transforms of the light curve shape moments 
Hl-fJ.3 show modulations in the transit shape with a period of 
3F or t,. The peaks around 0.189 1/day frequency (= l/P or (,/3) 
have 4-cr confidence. 



The semi amplitude of this variation is 21 ppm, well above 
the average noise level of 1 ppm in the Fourier spectrum. 
Harmonics of this frequency can also be detected up to the 
fourth order. This period has been found by Shporer et al. 
2011 and Mazeh et al. 2011 independently. The binned phase 
diagram of the residuals with 25.43 hour period is plotted 
in the bottom panel of Fig. [T] 

Both Mazeh et al (2011) and Shporer et al (2011) 
suggested a probable pulsational origin. Contrary to their 
proposition, we interpret it as the rotational period of KOI- 
13 A. The arguments for a rotational origin are the following: 

• The period is perfectly compatible with the expected 
rotational rate. Barnes et al. 2011 predicted a 22-22.5 hour 
rotational period for KOI-13 A, depending on the stellar 
mass. They did not estimate the error, but it is easy to 
calculate that errors in v sin i and in the stellar inclination 
results in ±3.9 hours. This prediction is in good agreement 
with the rotational origin of the new period. 

• Balona (2011) detected signs of stellar rotation in 20% 
of all A-F stars in the Kepler field caused by the granula- 
tion noise, and 8% of them also exhibit starspot-like features 
in the light curves. Typically, these stars exhibit a domi- 
nant peak with 10-100 ppm amplitude at periods less than 
3 days, the median is around 1 day for stars in the range of 
7500-10,000 K. Another diagnostic is the scatter level at fre- 
quencies below < 50/ day, which exponentially increases by a 
factor of ssl.6 toward low frequencies/high periods. The pe- 
riodogram of KOI-13 looks exactly as described by Balona 
(2011), and the folded light curve (Fig. 1, lower panel) is 
also compatible with a rotating A-type star with magnetic 
activity. 

• Stars similar to KOI-13 A can exhibit p or g mode pul- 
sations (5 Set and 7 Dor) or both (Uytterhoeven et al. 2011). 
However, the observed pulsations have periods less than one 
day, and several modes are observed with similar amplitudes. 
The general appearance of the frequency spectrum of KOI- 
13 does not resemble a pulsating star. 
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• Mazeh et al 2011 detected four harmonics of the 25.4- 
hour frequency in Q3 data, which we confirm here. The pres- 
ence of such harmonics is typical for stellar rotation (Balona 
et al. 2011). 

It should be noted that both we and Mazeh et al. 2011 
detected other frequencies between 1.5-2 c/day which are 
separated about equidistantly and have amplitudes of 4-7 
ppm. Their harmonics do not appear in Q2+Q3 data, and 
these peaks do look compatible with pulsation. 

The source of the 25.4-hour period is the host star, KOI- 
13 A. This is seen from the systematic modulation of the 
transit shapes. Because this period is in 3:5 ratio with the 
orbital period, every third transit occurs in front of the same 
stellar surface, and a modulation of the light curve shape 
is expected with a period of 3 transits. Since the individual 
light curves are too noisy for a direct comparison, we have to 
combine many data points and analyse their moments. Let 
us define fj, n , the n-th light curve moment of each individual 
transit as 

i£ {transit} 

where ti and fi are the times and occulted fluxes belonging 
to each data points, D is the transit duration and d is the 
calculated mid-time of the transit, based on the ephemeris 
in Borucki 2011. Once the moments are assigned to the in- 
dividual transits, a time series of moments can be analyzed 
in the standard fashion. In Fig. [2] we plot the periodograms 
of the first three lightcurve moments (note that the ampli- 
tudes are in relative units). The periodograms show a de- 
tection with 5.27 day period, with a significance of 3 — 4-cr 
for each light curve moments, confirming that the 25.4-hour 
period modulates the shape of the transit with a period of 

3Porb — 5-/~Vot • 
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Figure 3. Variations in the transit duration (top panel), the 
reciprocal transit duration (middle panel) and in the square of the 
impact parameter derived from transit shape (b 2 , lowest panel), as 
the function of the cycle number. The plot of the reciprocal transit 
duration is superimposed by the model of the best-fit linear trend 
together with the 3<r errors of the linear regression. 



2.1 The possible 5:3 spin— orbit resonance of 
KOI-13 A 

The 25.4-hour period is very close (within 0.1%) to the 5:3 
spin-orbit resonance with KOI-13. 01. Because of the large 
mutual inclination of 59° (Barnes et al. 2011), the longitude 
of the companion varies with changing velocity. Interest- 
ingly, the sinodical longitude of KOI-13. 01 remains exactly 
constant (within 1% fluctuations) on « 1/8 orbital arc sur- 
rounding the positions when the companion is at extremely 
high/low latitudes. Thus, KOI-13. 01 and KOI-13 A move as 
if they were in exact 1:1 resonance for 3 hours. 

To date, theories of spin-orbit resonance cover the fol- 
lowing cases: (i) spatial approximation with rigid bodies (e.g. 
Makarov 2011), where the body in resonance has little mass; 
(ii): resonant orbits of massive stars (e.g. Witte & Savonije 
2011); (iii) resonances in systems of compact bodies (e.g. 
Schnittman 2004). The case of resonance between the stellar 
spin and the orbit of a substellar companion is yet to be ex- 
plored theoretically but, if confirmed, it will be a significant 
finding in relation to the evolution of planetary systems. 



3 TRANSIT DURATION VARIATIONS 

The analysis of the light curve moments suggested a long- 
period time-dependent variation in the transit curve shapes. 
By fitting each transits individually, we concluded that the 
duration of the transit is gradually increasing (Fig. top 
panel). The significance of this finding is 3.85ct, with a false 
alarm probability of 0.012%, based on an MCMC estimate. 
The 3-er confidence interval of the fitted linear regression is 
plotted in the middle panel of Fig. [3] onto the distribution 
of the reciprocal transit length (which is the actual output 
of the light curve fit algorithm). We now describe the details 
of how the light curves were analyzed. 

First, every tr ansit was fitted i ndepe ndently using the 
analytic models of iMandel fc Agoll (120021 ). This fit invokes 
symmetric templates which cannot fit the known asymmetry 
of the light curves; but since the asymmetry appears simi- 
larly in each transits, we expect no time-dependent bias in 
the results. The free model parameters are: the transit time, 
the relative radius of the planet, the reciprocal of the half 
duration, Q/R*, and the square of the impact parameter, b 2 . 
These parameters are almost uncorrelated to each other (Pal 
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2008). Indeed, our analysis shows that the inverse duration 
C/R* follows a secular trend (see Fig. 3). 

The derived values of these parameters are displayed in 
Fig. [3] Since KOF13 did not exhibit a transit timing vari- 
ation, we can rule out short-term secular variations in the 
orbital semimajor axis of the transiting companion. There- 
fore, we conclude that the decrease in (,/R*, indicating a 
lengthening of the transit duration, is due to the increasing 
orbital inclination (i.e. due to decreasing impact parame- 
ter). The observed linear trend in C,/R* is d(£/i?*)/dt = 
(-31.6 ± 8.2) • 10~ 5 d" 1 cycle" 1 = (-17.9 ± 4.6) • 10" 5 d" 2 . 
Since the relation between (,/R*, a/_R* and b 2 is 



a } _ VI - b 2 ( C, 



KR.J> (2) 
(see also Ipil [20081 '). we can compute the time derivative of 
b as 



db 

dt 



1-b 2 



1 A [A. 

dt KR* 



(3) 



This was obtained by re-ordering equation @ and assum- 
ing that a/R* is constant. By substituting our best-fit 
b = 0.253 ± 0.020 into this equation, we had b = (-4.4 ± 
1.2) x 10" 5 d" 1 = (-0.016±0.004) y _1 . This is a tiny change 
in the transit duration, an order of magnitude smaller than 
the precision with which b can be determined from the light 
curve shape in the fitting procedure (Fig. 3, bottom panel). 
More precisely, the planet appeared to moved by only ~15% 
of its radius during the observed «100 transits, having no 
detectable effect on the light curve shape. Consequently, we 
reasonably neglected the oblate shape of the rotating star 
when converting transit duration to b. The stellar oblate- 
ness of a 2 Mq star with radius 1.7 Rq is 3% (Murray & 
Dermott, 1999), and less in an inclined projection. Therefore 
the bias due to stellar oblateness is 1-2% in the determined 
value of b, much below the accuracy of its determined value. 
In Section 4, we will show that b is compatible with the sec- 
ular perturbations caused by the oblateness of KOI-13 A, 
the rotating host star. 



4 INTERPRETATION 

Assuming M v = 9.2 Mj (see also Barnes et al. 2011), the an- 
gular momentum in the star and the companion are similar. 
Thus, the axes of both the planet orbit and the stellar spin 
are precessing around the total angular momentum axis, op- 
posing each other. To date, there is no exact theory for this 
case. In the observed examples, angular momentum of the 
orbit (double stars) or stellar spin (e.g. Mercury or lunisolar 
precession) dominates the other. One can study the preces- 
sion of the orbital plane and the stellar spin under slightly 
different assumptions. We will show here that these lead to 
compatible results, and give a satisfactory estimate of the 
stellar oblateness. 



4.1 Secular J 2 perturbations 

As known from the theory of satellite motions (|Kaulalll966h . 
higher order moments of the gravitational potential of a host 
body yield periodic and secular perturbations in the orbits 
of nearby companions. The external gravitational potential 
of an extended body can be expressed as 



n=2 



-J P n {cos8) 



(4) 



where M is the total mass, R is the equatorial radius, J n 
are constants and V n are the Legendre polynomials. The 
most prominent perturbation is caused by J2, due to the 
oblateness of the host body. MacCullagh's Theorem allows 
us compute J2 using 



®xx 



MR 2 



(5) 



where Q xx — Q yy ^ Q zz are the principal moments of iner- 
tia. It is known that a non-zero J2 results in secular pertur- 
bations in the angular orbital elements. Namely, the secular 
term in fi (argument of ascending node) is computed as 



dfi 
~dt 



\J2r1 



(-Y 2 - 

\RJ (1 



cos tp 
(1 - e 2 ) 2 



(6) 



Here n denotes the orbital mean motion, a is the semi-major 
axis and e is the orbital eccentricity. 

It is known from vector geometry, that if an unit vec- 
tor n precesses around the unit vector p with an angular 
frequency of ujq, the time derivative of n will be 



h = ojo(p x n) 



In our case, n is the unit vector parallel to the orbital angular 
momentum of the transiting body. Since the components of 
n are 



sin i cos Q 
sin i sin Q 
cos i 



(8) 



and the only observable quantity is n z = cosi, we can write 

h z = u) (p x n y — p v n x ). (9) 

By substituting p x = sin i p cos fi p and p y — sin i p sin Sl p and 
the components n x and n y from equation @, we obtain 



dcosi 
dt 



= u>o sin i sin i p sin A, 



(10) 



where fi p is the ascending node of the stellar equator, and 
A is the longitude of the planet's ascending node, relative to 
that of the ascending node of the stellar equator by defini- 
tion. 



4.2 The inferred oblateness of the host star 

Assuming a circular orbit for the transiting companion and 
substituting the above relation and equation (??) into equa- 
tion ([6]), and by taking uq = dtt/dt as the precession rate 
inducted by J2, we finally obtain 



dcosi 



df 



-(E) >< 
x (cos i cos i p + sin i sin i p cos AS!) x 

x sin i sin i p sin AQ. 

Since transits are observed, we can say here that cos i 
sini « 1. In addition, b = (a/i?*) cosi, thus the above equa- 
tion can be rearranged to give b as 



-Jin 



(12) 
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For J 2 we obtain 

J 2 sin 2 i p sin A cos A = (3.8 ± 1.0) • 1CT 5 . (13) 

Substituting the derived stellar parameters we find J2 = 
(2.1 ± 0.6) x 10~ 4 , and dtt/dt = (3.4 ± 0.9) x 10~ 5 / da y- 

Is this result compatible with a model star with the 
same mass, radius and rotational rate as KOI-13 A? To 
check this, we calculated the moment of interia of n — 3 
polytrope models that were distorted uniform stellar ro- 
tation. The distortion of a rotating star is (R — r)/R = 
Q 2 R 3 (2GM)- 1 + 3J 2 /2 ||Stbdl2004h . where R and r are 
the radii at the equator and the poles, respectively. The 
3i term we derived for KOI-13 is two orders of magnitude 
less than the rotational term and can be neglected. The 
local distorions were calculated at at each internal point 
and then a ^-contraction was applied. We assumed dz/z = 
-Q. 2 r 3 (2GM r )~ 1 = -3D. 2 /(8nG(p)r), where M r and (p) r 
are the mass and the mean density inside the radius r. This 
model is a first approximation that neglects higher-order dis- 
tortions, but is enough precise for a rough estimate because 
the distortions are everywhere smaller than 3%. The result- 
ing values are Q zz = 0.07760 MR 2 , @ xx = 0.07743 MR 2 , 
J2 = 1.7 x 10 -4 . Thus, the expected value of J2 is in the 
range suggested by the theory of secular perturbations, sup- 
porting our interpretation. 



ular variations can be detected more e asily in TDV than i n 
transit timings for short time coverage (|Pal fc Kocsisll200gft . 
An argument against the reality of TDV can be that we 
fitted symmetrical templates to a transit with known slight 
asymmetry. However, we consider that the detected trend is 
not an artifact of this kind, because (i) there is no visible 
time-dependent trend in the asymmetric part of the light 
curves; (ii) the errors contain the ambiguity introduced by 
the asymmetric part, while the detection is still significant, 
and (iii) it is rather unlikely that the varying asymmetries 
affect only the transit duration, and leave the other param- 
eters unchanged. 

In interpreting the cause of transit duration variations, 
our hypothesis of secular J2 perturbations is a probable sce- 
nario, but other processes cannot be excluded (e.g. pertur- 
bations from an unreported outer planet; the presence of a 
moon around KOI-13. 01). However, the stellar rotation is 
the most plausible reason for TDV, because it does not in- 
voke other bodies into the explanation, and the inferred J2 
is fully compatible with the expected value. In other words, 
one would predict secular J2 perturbations of KOI-13. 01 at 
the rate detected, taking into account the rapid rotation of 
the host star and the close proximity of its substellar com- 
panion. 



4.3 Precession rate 

Applying the framework developed for lunisolar precession, 
the precession rate of the star can be formulated as 

d£l p 3GM P (Q ZZ — 0„) cos ip 3 GMpCosip . . 

~dt = 2a 3 Q zz n = 2 a 3 nf3. ( ' 

Here M v is the mass of the planet and /3 = 0.07760 for KOI- 
13 A, as we derived in the previous section. Evaluating Eq 
14, and applying the value of J2 = 2.1 x 10~ 4 we find that 
dQ. p /dt is in the range (3.81-1.67) x 10 _J /day, for a compan- 
ion of mass between 4-9.2 Mj. Therefore, dQ. p /dt has the 
same order of magnitude that we derived for d^l/dt. The pic- 
ture of a planet's plane precessing with « 500 year period, 
and a star's orbital axis precessing with the same period, is 
self-consistent for KOI-13. 



5 SUMMARY 

With high-precision Kepler photometry, we detected photo- 
metric and dynamical effects of stellar rotation in the ex- 
oplanet host KOI-13. From photometry alone, we have de- 
tected the following three phenomena for the first time: 

• Stellar rotation in a probable resonance with the orbital 
period of the close-in substellar companion; 

• Long-period transit duration variations; and 

• Precession of the orbital plane of an exoplanet candi- 
date. 

Variations in transit duration were found from the fit of 
the consecutive light curves to transit shape models. Simul- 
taneously, we fitted other model parameters such as the rela- 
tive radius of the planet, the impact parameter and the tim- 
ings of transits, and no variations were found in these quan- 
tities. The negative TTV detection is plausible because sec- 
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